#' @param times (`numeric()`)\cr
#'  If `integrate == TRUE` then a vector of time-points over which to integrate the score.
#'  If `integrate == FALSE` then a single time point at which to return the score.
